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SUMMARY 

A correspondence principle has been developed for treating the steady- 
state propagation of waves from sources moving along a plane surface or inter- 
face. This new principle allows one to obtain, in a unified manner, explicit 
solutions for any source velocity. To illustrate the correspondence principle 
in a particular case, the problem of a load moving at an arbitrary constant 
velocity along the surface of an elastic half-space is considered. 

INTRODUCTION 

Certain problems in the linear theory of wave propagation are of fundamen- 
tal importance to a wide variety of fields. One of these is the response of a 
plane interface between two different materials to moving transient sources of 
disturbance. The reflection and refraction of plane transient waves at an in- 
terface (refs. 1-3) and the generation of waves from specified sources moving 
at a constant velocity along an interface (refs. 4-7) are two important exam- 
ples of this type of problem. In such cases it is usually assumed that the 
surrounding media is in plane motion and that a steady-state wave pattern ex- 
ists relative to an observer moving with the source of disturbance. The re- 
sulting two-dimensional steady-state boundary value problem can then be solved 
either by transform techniques or by the use of complex function theory and 
the method of characteristics (see refs. 8 and 9) . 

In many problems, however, both of these traditional methods are very in- 
efficient. This is because it is necessary to pose and solve separately the 
special cases when the source velocity is less than or greater than each of the 
characteristic wavespeeds in the surrounding media. This paper demonstrates 
that it is possible to treat all such special cases in a simple, unified manner 
through the application of a newly developed correspondence principle. In ad- 
dition, this correspondence principle leads to a new and direct representation 
for the general solution of steady-state interface problems. 

PROBLEM STATEMENT 

Consider two homogeneous, isotropic se m i-infinite media (either fluid or 
solid) which are in contact along the plane y = 0 and which contain disturbances 
traveling at a constant velocity U in the negative x-direction. We assume that 
these disturbances are uniform in the z-direction and that a steady-state motion 
exists in the semi-infinite media. Under these conditions, the governing equa- 
tions of motion in the two media reduce to (ref. 8) : 
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(m - 1,. ., j) 


(1 - U 2 /c 2 * )3 2 <J> /3x 2 + 3 2 <f> /3y 2 = 0 
mm m 


(i) 


in a set of moving coordinates (x,y ,z) defined by x = x + Ut, y = y, z = z. In 
equation (1), (|> m = <j> m (x,y) and c m are the displacement potentials and their 
corresponding wavespeeds, respectively, for the two media. 


Along the surface y = 0, the <f> m must satisfy a certain set of boundary or 
continuity conditions which, in general, can be written in terms of the second 
order partial derivatives of the <}> m as 

^^h/ 9 * 2 ’ 9 \/ 9x9y ’ ^n/ 9 ^ = P n^ (n = 1 , . . . j ) (2) 

where the E n also depend on the source speed U and the material properties of 
the two media and are linear functions of their arguments. The vector - 
(P n (x) } is determined by the values of the source disturbances along the plane 
y = 0 and is assumed to be given. 


Since the linear operators which appear in equation (1) are hyperbolic if 
U > c m and elliptic if U < c^, the steady-state solutions of these equations 
will depend on the relative size of U and the wavespeeds c^ Consider first 
the "totally supersonic" (TSS) case (i.e. where U > Cju is satisfied for all m 
in equation (1)). 

Totally Supersonic Case 

In the TSS case, the general solutions to the equations of motion (1) can 
be written as 

'6 = F (x - 6 | y | ) (3) 

mm m 

where g m = (U^/c^ - 1)1/2 an< j the bars denote "absolute value of". In equation 
(3), solutions of the type F m (x + 3 m |y|) have been rejected since they repre- 
sent disturbances traveling in the positive x-direction and, hence, would vio- 
late the "radiation conditions" (ref. 8) . 


The second order partial derivatives of the <j> m then become 

p 9 " 

94 /3x Z = F 
m m 

3 2 4 /3x3y = - sgn(y)6 F 

m mm 


(4) 


2 2 
3 * /3y 
r m J 


= g 2 F 
m m 


ii 


where F m denote the second derivatives of F m with respect to their arguments 
and sgn stands for "sign of". Placing equation (4) into equation (2) thus 
yields a set of linear equations in the F m on y = 0, which using the summation 
convention can be written as 
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(n — lj • • • 3 ) 


it 


A F 
nm m 


(x) = P n (x) 


(5) 


Here, the A is a real jxj matrix involving only the material properties of the 
given media and the disturbance speed U. Assuming that the matrix is non- 
singular, we can then solve for the F m , obtaining formally 


F <*) = 


m 


A -1 P (x) 
ran n 


( 6 ) 


and these results can be continued into the two media, using equation (3), to 
give: 


m 


" -i 

= A P (x 
mn n 



(7) 


The stresses in each media can be obtained directly from equation (7) since 
they are simply linear functions of the <|> m . The displacements or velocities, 
however, must be found from these results by a single integration. This TSS 
case is of fundamental importance for steady-state interface problems of the 
type we have been discussing because it contains implicitly, through equation 
(6), the solution for all other cases when this equation is interpreted in an 
operational sense. To prove this we now consider the general case. 


General Case 


Consider the general case when U/c m < 1 for m = 1,... k, where k j . 
Then the governing equations (1) are hyperbolic for m > k and elliptic for 
m £ k, and their general solutions are 


<f> m = Re{G m (x + igjy|)} m < k 

<f> = F (x - gjy|) k < m < j (8) 

m m m ” 

where g m = (1 - U 2 /^) 1 / 2 , Re denotes "real part of", and the G m for m <_ k are 
analytic functions of the complex variables x + ig m |y|. For m > k, the second 
order partial derivatives of the cf> m are again given by equation (4) , For m _< k, 
we now obtain instead 


2 2 " 

3 4> /3x = Re{G } 

m m 

9 2 i|> m /3x3y = - sgn(y) g Im{G } 
m m m 

9 <f> /9y = - gf Re(G } 

m mm 


(9) 


where Im denotes "imaginary part of". However, on the boundary y - 0 the real 
and imaginary parts of these G m satisfy a pair of Hilbert transforms 

Re{G } = H[Im{G "}] 
m m 

Im{G } = - H[Re{G "}] 
m m 
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where the Hilbert transform is 


H[f] = 1 /tt f fd?/(5 - x) 

—00 

and the integral is understood to be taken in the principal value^sense. Hence, 
the partial derivatives on y = 0 can be written in terms of Re{G m } only as 


9 2 <j> / 9x 2 = Re{G "} 
m m 

8 2 cj) /9x9y = sgn(0+)g H[Re{G }] 
m. — m m 

2 2 - 2 ft 
3 4. / 3y Z = - g Re{G } 
m m m 


( 10 ) 


Note now that on the boundary y = 0, equations (4) and (10) will be identical 
if we make the following replacements in the TSS expressions (4) for m k. 

g -> -ig (11) 

mm 

and identify the F^" in the TSS case with the G m " in the general case through 
the additional replacements for m £ k given by: 


If If 

F Re{G } 
m m 


iF ^ H[Re{G ' }] 
m m 


( 12 ) 


Thus, on the boundary y = 0 there is a one-to-one correspondence between the 
complex-valued TSS problem obtained by making the substitutions given by equa- 
tion (11) and the general case problem if, as equation (12) shows, the appear- 
ance of the imaginary number i in the TSS problem is interpreted as represent- 
ing the Hilbert transform operator in the general case. This correspondence 
also means that the complex-valued matrices A and A“^, which result in equations 
(5) and (6) from the substitutions given by equation (11), must be interpreted 
as representing matrix operators in the general case. In particular, breaking 
A~1 into its real and imaginary parts, we have 


A ^ = a + ib ■* a +b H [ • ] 
mn mn mn mn ran 


(13) 


where a mn and b mn are both real. Using this result and equations (6) and (12) , 
we see that on y = 0 the general case solution is given by 


F (x) = aP (x) + b H[P(x)] k < m < j 
m mn n mn n — 


Re(G (x)} = a P (x) + b H[P (x)] m < k 
m - mn n mn n — 


(14) 


Since the general solutions for m > k are constant along the real characteris- 
tics x - g m |y| (see equation (8)), the F m " can be continued directly into the 
adjacent media and the general case solution written as 
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+ " = a P ( x “ | y | ) + b H[P (x - 6 |y|)] k<m < j (15) 

m mn n n mn n rx — 

For m _< k, however, our problem consists of finding the functions G m "(x + i-8 m |y|) 
which are analytic in the upper half of the complex-plane and whose real part is 
given on the real axis by equation (14) . This is a standard problem in analytic 
function theory whose solution may be written as 

°° -1 

G = 1/iir f AjT (£)d£/(£ - z ) m < k (16) 

m mn n in — 

— oo 

provided the integral converges, where z = x + iS m |y|. As before, the stresses 
can be obtained directly from equations (15) and (16) although a further inte- 
gration is necessary for displacements and velocities. 

With the general solutions given by equations (15) and (16), it is now 
particularly easy to obtain the solution to steady-state interface problems for 
arbitrary source velocity. All that is needed is the inverse matrix A“1 from 
the TSS case solution. In the general case this matrix becomes complex-valued 
when the substitution in equation (11) is made. A simple algebraic decomposition 
of A~1 into its real and imaginary parts for each special case then gives the 
necessary matrices for the expressions in equations (15) and (16). To illustrate 
the use of this method we now consider a particular problem, 

MOVING LOAD ON A HALF-SPACE 


A number of authors (refs. 4-7) have previously considered the response of 
an elastic half-space to loads traveling at a constant velocity on the plane 
surface. Here, we will solve for the waves generated in the half-space y >_ 0, 

-oo < x < °°, -°° < z < 00 by a moving distributed load of intensity P(x) in the mov- 
ing coordinates x=x+Ut, y = y, z=z (figure 1). Then the normal stress, 
tyy, and shearing stress, t . on the surface are given by 

t = -P(x)sin0 

yy 

t = P(x)cos9 
xy 

where 0 is the angle between the direction of the applied load and the half- 
space surface. In this case there are only two displacement potentials <{>-|_ an< l 
(f>2> which correspond to dilatational and shear wave disturbances, respectively, 
and two corresponding wavespeeds c^ and C 2 - Application of the boundary condi- 
tions (17) yields the matrix A and vector P^ given by (ref. 8) : 


(17) 


(M 2 - 2) -20 2 


— —i 

-Psin0/p 


P = 


-28 1 “(M^ - 2) 


Pcos8/p 


(18) 


where p is the shear modulus and M 2 = ^/ c 2 ' Then the inverse matrix A ^ is 
given by 
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(19) 


(M^ - 2)/D -28 2 /D 



|^-2 8 1 /D -(M 2 - 2)/Dj 

2 9 

where D = (M 2 - 2) z + 48^32 • Table 1 shows the breakdown of this inverse matrix 
into its real and imaginary parts. In that table = (M^ - 2)^ + 168x32 an <l 
D 2 = (M^ - 2)^ - 48^82* When those results are placed back into equations (15) 
and (16), the problem is then formally complete. To illustrate the use of these 
expressions for a particular loading, consider the case of a moving concentrated 
line load, i.e. P(x) = P<5(x) where P is a constant and S(x) is the Dirac delta 
function. Then we obtain: 


Totally Supersonic Case (U > c^ > c^) 

It 

<(> 1 = (-arsine + a 12 cds0)P6(x - 8^)/^ 

ft 

$2 = (-arsine + a 22 cos0)P6(x - 8 2 y)/y 


Transonic Case (c 2 < U < c^) 

(a^sin0 - a^ 2 cos0)P/iiryZ^ + (b^^sin0 - b^ 2 cos0)P/iryZ^ 
( -a 2 l s in0 + a 22 cos0)P6(x - 3 2 y)/h 
+ (b 21 sin0 - b 22 cos0)P/iry(x - 8 2 y) 

Subsonic Case (U < c 2 ) 

II 
II 


a^sin0 - ib^ 2 cos0)P/i7rjjZ^ 


-a 22 cos0 + ib 2 ^sin0)P/iirpZ 2 


! I 



Similar results to these have been derived by the traditional complex variable 
and characteristics approach in the treatise by Eringen and Suhubi (ref. 9). 


CONCLUDING REMARKS 


The correspondence principle developed above has led to a new unified form 
of the solution for steady-state interface problems (equations 15 and 16) which 
can be efficiently used to treat a number of problems. In addition, this prin- 
ciple clearly demonstrates the close relationship that exists between the struc- 
ture of the general solution and the TSS case. This relationship is currently 
being extended to steady-state problems in anisotropic media. 
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